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Abstract 

Strongly consistent and asymptotically normal estimators of the Hurst parameter of solutions of stochastic differential 
equations are proposed. The estimators are based on discrete observations of the underlying processes. 
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1 Introduction 

Recently long range dependence (LRD) became one of the most researched phenomena in statistics. It appears in 
various applied fields and inspires new models to account for it. Stochastic differential equations (SDEs) are widely 
used to model continuous time processes. Within this framework, LRD is frequently modeled with the help of SDEs 
driven by a fractional Brownian motion (fBm). It is well known that the latter Gaussian process is governed by a 
single parameter H £ (0,1) (called the Hurst index) and that values of H in (1/2,1) correspond to LRD models. In 
applications, the estimation of H is a fundamental problem. Its solution depends on the theoretical structure of a model 
under consideration. Therefore, particular models usually deserve separate analysis. In this paper, we concentrate on 
the estimation of H under the assumption that an observable continuous time process (A't)tg[o,T] satisfies SDE 

Xt=i+ f f{Xs) ds + r giXs) dBf, t £ [0, T], (1.1) 

Jo Jo 

where T > 0 is fixed, / is an initial r.v., / and g are continuous functions satisfying some regularity conditions, 
and {B^)t^[o,T] is a fBm with the Hurst index 1/2 < H < 1. Our goal is to construct a strongly consistent and 
asymptotically normal estimator of the H from discrete observations Xtj,..., of trajectory Xt, t £ [0, T], 

We consider two cases. First, we examine the case when g is completely specified. Next, we relax this restriction 
and allow g to be unknown. Such situation may appear when g depends on additional nuisance parameters. In both 
cases, boundedness of 1/g plays an important role and is assumed to hold. 

To our best knowledge, so far only several studies investigated this question. The pioneering work has been done 
by [3] of Berzin and Leon as well as lecture notes [4] with references therein. [12], [14] and [15] were also devoted to 
the problems of the same nature. However, all of these works focused on the strong consistency. The present paper is 
a generalization of [13] where a special case of (1.1) was considered. 

The paper is organized in the following way. In Section 2 we present the main results of the paper. Section 3 is 
devoted to several results needed for the proofs. Sections 4-5 contain the proofs of the main results. Finally, in Section 
6 two examples are given in order to illustrate the obtained results. 


2 Main results 

To avoid cumbersome expressions, we introduce symbols O^,, Oi^. Let (Vi) be a sequence of r.v., ? is an a.s. non-negative 
r.v. and (a„) C (0, oo) vanishes. Yn = Ooj(a„) means that |y„| ^ c • a-n', Yn = means that \Yn\ ^ ■ bn with 

bn = o{a„). In particular, Yn = Ou,{l) corresponds to the sequence (Yn) which tends to 0 a.s. as n ^ oo. 

‘Corresponding author. E-mail: kestutis.kubilius@mii.vu.lt 

tThis research was funded by a grant (No. MIP-048/2014) from the Research Council of Lithuania. 

“E-mail: viktor.skorniakovamif .vu.lt 


1 



Let TTn = {t^ , k = 0,..., in}, n > 1, N 9 in t oo, be a sequence of partitions of the interval [0, T], If partition vTn 
is uniform then ^ for all k G {0,. .., in}- If in = n, we write instead of rj}. In order to formulate our main 

results, we state two hypotheses: 


(H) 

(Hi) 


AX^n = Xnn - = 0„ (d"-") , k = 1, . 

= X^n - 2Xn" + = alX^n 

^k — 1 ^k — 2 ^ ^ ~k — 




Ouj {d 


2(H-e) 

n 




k — 2,..., in, 


( 2 . 1 ) 


for all e € {0, H — i), where dn = maxi^fc^i„ {tJ} — 

Theorem 2.1. Assume that solution of Eq. (1.1) satisfies hypothesis (Hi). Moreover, let the funetion g is known, 
Lipschitz-continuous and there exists a random variable c such that P(<; < oo) = 1 and 


I M 

te[o,T] \9\.^t)\ 


^ ? a.s. 


( 2 . 2 ) 


Then 


H, 


( 1 ) 


H a.s., 


2^ In - -H)^ Nifi- a%) for H G (1/2, 1), 


where a% is a known variance defined in Subseetion 3.1, 


=¥>“V 






for n > T, 


Tn,T{x) =(^~j (4 — 2 ) and (p„ j, denotes an inverse of ipn,T, x G (0,1), n > T. 


□ 


It is natural to try to drop restriction of the known g. For this purpose we need several additional definitions. 
Assume that the process X is observed at time points ■:^T, i = 1,..., mn, where mn = nkn, and kn grows faster than 
nlnn, but the growth does not exceed polynomial, e.g. fen = nln® n, 6 > 1, or kn = n^. 

Denote 

U U 

^ 2 


IFn.fc = 


J = -fcn + 2 

where I ^ k ^ n — I and sf = -i—T. 

J rrin 




l = - fc „+2 


Xs"+e — 2Xs>» +e + Xs" „+e 

3 k j—l'k J—2 k 


)■ 


Theorem 2.2. Assume that solution of Eq. (1.1) satisfies hypotheses (H) and (Hi). Moreover, let the funetion g is 
Lipsehitz-continuous and there exists a random variable c such that P(c < oo) = 1 and inequality (2.2) holds. Then 


—> H a.s., 

^ >-Y-\ d 


2^ In - (#(2) X(0; o\) for H G (1/2, 1), 


where 


H, 


( 2 ) _ 


+ 


1 


2 In k, 




and -I is a known varianee defined in Subsection 3.1. 


3 Preliminaries 


3.1 Several results on fBm 

Recall that fBm with the Hurst index H G (0,1) is a real-valued continuous centered Gaussian process with 

covariance given by 

tjv-dRdRx 1 / 2ir , ,2H I, |2ir\ 

E{Bt ^ -[s +t -|t-s| ). 

For consideration of strong consistency and asymptotic normality of the given estimators we need several facts regarding 


Limit results. Let 


Vn.T = 


r2R(4_ 22R) 


E 


H\2 
tu) ’ 




1 

2 ■ 


2 



Then (see [2], [10], [1]) 


V„,T ->■ 1 a.s., 

n—>-oo 

where 

Holder-continuity of . It is known that almost all sample paths of an £Bm are locally Holder of order strictly 
less than H, H £ (0,1). To be more precise, for all 0 < e < H and T > 0 there exists a nonnegative random variable 
Ge,T such that E(|Ge,T|^) < cxD for all p ^ 1, and 

|Bf-Bfl ^Ge,T|i-s|"-" a.s. (3.1) 


for all s, t € [0, Tj. 

3.2 Concentration ineqnality 


Let 


^k.n — 


T»V4- 22» 


^(2)bH 


t? > ^ Ey,>yfc,„, j,k = 2,...,n. 


Note that = pnij — k). In the sequel we make use of the following modified version of an inequality of concentration 
from [6]. 

Lemma 3.1. For all z > 0 and any H £ (0,1), 

P( >2^ ^ 2 exp f - ^ 

Proof. Let x = sup^g^p jj Sjgz \pH{j)\- Following an argument of the paper [6], one gets bound 

n n 

2 (n-i)-i ^ 2 (n-i)-i in,„|-|y,>|-|p«(i-fc)| 

k,j=2 kyj=2 

n 

«:2(n-l)-i Yl„\pHij-k)\ 

k,j=2 

n n 

^2{n - = 2(n - 1)“^;^ {yI„ - l) + 2>r 

k=2 k=2 


Thus (see 


In paper [5] it was proved that 


and 




> 2 : < 2 exp < — 






lez 


for H^l/2, 


! + ■ 


= 2, for H ^ 1/2, 


>< = sup E \Pn{j)\ = lim E IPw(i)l = - . 


This yields the required inequality. 
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4 Proof of Theorem 2.1 

Proof of Theorem 2.1. Observe first that 

is continuous and strictly decreasing for n > T. Thus, it has an inverse for n > T. By hypothesis (Hi), 


(l)^ 


Tn,T{H, 
Tn,T{H) 




y2ff(4_ 22ff) 

IH — \ Tl . . 

>211(4 _ 22ff) ^ ( n«-3" ) "" + O- ( ^H-3e ) 


1 


for 3s < H. Therefore 




1 a.s. as n —>■ oo. 


Tn,T{H) 

Using the same arguments as in [13] it is possible to prove that the estimator is strongly consistent and asymp¬ 
totically normal. 

5 Proof of the main Theorem 

Before presenting the proof of this theorem, we give two auxiliary lemmas. 

Lemma 5.1. Let 

2H k 


Vn,T{k) = 




2fc„r2ii(4_22ii) 




j = - fcn -|-2 


The following relation holds: 


max \Vn,T{k) - l\ = Ocu, , 

' \\kn 


Proof. By self similarity and stationarity of increments of fBm, 

kn 


Vn,T{k) = 


o2H-l_2H 
d 2 


22'»fc„(4-22«) 


l=-fcn+2 


^ BgiL - + B 


H 

.? + fcr 7.-2 


(2m„)^ 




(2fc„) = 




2211(4 _ 2211) V 

3 —2 


d {2kn 




H 

j-2 


+ 1 

2 


4_22h 


J=2 


^ B^-2i3f-i+B> =U2fc„.i. 


Therefore, 

P ^|Un,T(fc) — l| > 5 )<Ep(I> T(fc) - l| > <5) nP(|U 2 fc„.i - 1| > 5) for all 5 > 0. 

Put Vn,T = -^^Vn.T- Note that 

|V2fc„,l — 1> |V2fc„,l — l| -I- —— . 

2kn 


4 





Let ((5n) be a sequence of positive numbers such that 4' 0 n ^ oo and < ^n- By Lemma 3.1, 


P (|V2A;^,1 — 1| > 2^n) ^ P ^|P2A;„,1 “ 1| + ^ 2(5nJ ^ P (^\V2kn,l “ l| > 

= P (V2fc„-iy2fc„,i - 1| > 5„V2k„ - l) < 2 exp I 
Set 6n = y^fl 2fcn-i ~’ k„ ^ nlnn, then 


f(<5n + l) 


P (|P 2 fc„.i - 1| > 25„) ^ 2 exp < — 


3a In n 


32 >/a 


In n 


+ 1 


2n In n — 1 

If a ^ 3 and n ^ 2, then 5„ > k^^. Moreover, P (|P 2 fc„,i — 1| > 2(5„) 2n“® for a > 32 and n large enongh. Therefore 

series P(niaxi^fc^n-i I Vn,T(fc) — l| > 5n) converges and by the Borel-Cantelli lema maxi^fc^n-i I Vn,T(fc) — l| - > 0 

n n->oo 

a.s. 

Lemma 5.2. Assume that function g is Lipschitz-continuous. If e < {H — l/2)/3, then for each fc = 1,... ,n — 1 

n2H / A o2H\ 


W„,k g (Xtj) (2fc„)2»-i rnlff j + 

Proof. Step 1. By hypothesis (Hi), 


kn 




, j =-k^ + 2 ,... ,k„. 

\ TTln / 


Next, note that 


\g\Xt) - g\X,)\ =\g{Xt) - g{X,)\\g{Xt) + g{X,)\ ^ L\Xt - X,\\g{X,) + g{Xt)\ 

<:2L snp \g{X^)\\Xt - Xs\ 

uG[0,T] 

where t > s and L is Lipschitz constant. Thns, hypothesis (H) with a.s. continnity of t i—>■ g{Xt) lead to 

Step 2. Assume that e < {H — l/2)/3. Step 1 and Lemma 5.1 yield 


j = —fcn+2 j = — kn+2 


3(H-e) 

mn 


3 — ~fcn+2 j — — fcn+2 

kn 


+ Ou. 


3(H-e) 


m. 


_2,^ w 2fc„r""(4-2"") 


=9\X..) 


2fc„r2"(4- 2^") 


P„.T(fc) + 0^ 


2(H-e) 

n-" ^run 


+ 0.. 


\/fcn Im 


+ o^ 


+ 0. 

kr, 


3(H-e) 




Conseqnently the proof of lemma is completed. 
Proof of Theorem 2.2. Put 


. 2 ^ (A(^)A.n) 

nk^-^h 
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It follows from (2.1)-(2.2) and Lemma 5.2 that 


Sn.T =- 


nki 






5" 


2fcnT2^(4-22JJ) _|_ Q yfc;nS' 




+ 0. 


nfc„ 


I^tE 




fc=2 

v^2H 


HA.Z^I^ + 0^ (A^)+0^(. 

(A(2)B,f)" + 0. (n-3("-)) 


■E' 


nfc-T-(4 - 2-) ^,1 + 0^ (yi^) + O. {J^) 

I4.r + 0^(n-("-^")) 

^ (\/W) + (;r^) 

Term 0„ ^ ) vanishes provided e > 0 is small enough. Convergence Ki.t “'”'> 1 implies that Sn^r > 1. 

Hence 

^(2) ^ ^ ^ hi^ ^ 


2 In/Cn n—foo 

/'2') 

To prove asymptotic normality of the estimator Hn ' observe that 


Vn {Sn,T — l) =A 


K,t- 1 + 0.. JJp +0, 


+ Oi^ij 


-(H-Se)'- 


^ + o^{'JW)+o^{^) 

y/n{Vn,T — 1) 


o- (\/^) + o- + o- (» 


-(H-1/2-3£) ' 


+ 


for £ > 0 small enough. Now apply Slutsky’s theorem and limit results of Section 3.1. 


-^N(0,al,) 


6 Examples 

As mentioned previously, in this section we present two examples of applications of the obtained results. The hrst one 
deals with a general form of equation (1.1) and relies on certain restrictions on functions / and g. The second one 
describes a particular model which formally does not fit into the scope of the first one. 


6.1 Example 1 

In order to present an example, we need several facts on variation. To make the paper more self-contained and the 
structure clearer, the mentioned facts are briefly reminded in subsection 6.1.1. For details we refer the reader to [9]. 

6.1.1 Variation 

Fix p > 0 and —oo < a < b < oo. Let jr = {{a;o, • •., Xn} \ a = xo < ■ ■ ■ < Xn = b,n ^ 1} denotes a set of all possible 
partitions of [a, b]. For any / : [a, b] —>■ R define 

n 

Xpif-, [a, 6]) = sup^|/(a;fc) - f{xk-iW, Vp{f-, [a,b]) = [a,b]), 

^ k=i 

Wp{[a, b]) = {/ : [a, b] ^ R I i>p(/; [a, b]) < oo}, CWp{[a, b]) = {/ € Wp([a, b]) | / is continuous}. 

Recall that Vp is called p-variation of / on [a, b] and any / in >Vp([a, b]) is said to have bounded p-variation on [a, b]. 
For short we omit an interval [a, b] in the notations introduced above whenever there is no ambiguity. Below we list 
several facts used further on. 

• /!—>■ Vp{f) is a seminorm on Wp; Vp{f) = 0 if and only if / is a constant. 

• f €Wp^ sup,^e[a.6]l/(a;)l < oo- 
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• f,g fg e Wp. 

• q > p ^ 1 ^ Wp C Wq. 

• Let / € Wq, h € Wp with p,q G (0, oo) such that 1/p + 1/q > 1. Then an integral j*’ f dh exists as the 
Riemann-Stieltjes integral provided / and h have no common discontinuities. If the integral exists, the Love- 
Young inequality 


f dh - f{y)[h{b) - h{a)] 


^ Cp,qVq{f)Vp{h) 


holds for all y G [a, &], where Cp,q = C(P ^ + 9 and (/{s) = 5I]ri>i ^ Moreover, 


Vp( I fdh-,[a,b]] <Cp,qVq,^{f)Vp{h), 

^ a ' 

where Vq^^oif) = Vq{f) + sup^gj^ ^ l/(®)l- Also note that Vq^oo is a norm on Wq, q^ 1. 

Remark 6.1. The left-hand side of (3.1) can be replaced by [-Sif])) ^ {d,H — i), i.e. 

VHe{B^-, [s,t]) < Ge,T\t - s|^“^ a.s. (6.1) 

with Ge,T of (3.1). 


6.1.2 Assumptions and properties of solution of SDE 

Let (Rt^)tG[o,T], H G (1/2, 1), be a fixed £Bm defined on some probability space (fl, J-, P, F). Let a G (-^ — 1; l], 
^i+a(u) = R I /j' exists and sup^|/i'(a;)| -I- sup,,,^^ < oo}. Assume that / is Lipschitz and 

g G C^^“(R). In such case there exists a unique solution of (1.1) having the following properties: i) {Xt)t^[o,T] is F 
adapted and almost all sample paths are continuous; ii) Xq = ^ a.s.; iii) Vp > ■^, P{Vp{X-, [0,T]) < oo) = 1 (see [8], 
[16], [17] and [11]). 

Lemma 6.1. Let X satisfies (1-1), e G (0,77— ^). There exists a.s. finite r.v. L^.t such that 

VH,{X-,[s,t]) ^Le,T{t-sf^-‘^\ O^sKt^T. (6.2) 

Proof. Since Vp is seminorm and non-increasing function of p ^ 1, inequalities of subsection 6.1.1 give bound 

VHe{X-,[s,t]) \f{Xu)\du + Vp{^j^ g{Xu)dB//:-,[s,f\ 

^ sup \f{Xu)\{t-s) + GH^,HAvH,{goX.-,[s,t]) + sup |p(X„)|] [s,t]) 

tie[0,T] Y uS[0,T] J 

KLe,T{t—s)^ ^ 

with Le.T = sup„g[p_j,]|/(Xi.)|T^""^+® -I- Gh^.h^ (VH^{g o X.-[s,t]) -I-sup^gjo y]| 3 (X„)|^ Ge,T and Ge,T of (6.1). A. 
s. continuity of t i—>■ At together with continuity of / and g implie a.s. boundedness of sup^g[Q y]|/(Au)| and 
sup^gjQ y]|(;(A„)|. It is not difficult to show that I4r^ (g o A.; [s, t]) ^ sup,^|(j'(m)| (A; [s, t]). Hence, continuity 

of p'and a.s. boundedness of (A; [s, t]) guarantees a.s. boundedness of (p o A.; [s, t]) along with that of Lj^t. 

Lemma 6.2. Let X satisfies (1.1), e G {0, H — ^), and conditions stated above are true. Then the following relations 
hold: 

AA,~ =0,,(d"-"), fc = l,...,i„, (6.3) 

A(2>A,n =g(A,n_JA(2)B" + 0^{dT-^^), fc = 2,... ,i„, (6.4) 

where d„ = maxi^fcs;i„ {xk — 

Proof. Let a sample path t i—>■ At be continuous. We first prove (6.3). Note that 

f{Xs)ds + r [p(A.) - fl(A,n_ J]dBf + p(A,n_ JAB/^n. 
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An application of inequalities (6.1)-(6.2) together with continuity of /, g,g' and mean value theorem yield 


\f{Xs)\dsS^ sup \f{Xu)\{Tk - Tk-l) = Ou:{dn), 
u6[0,T] 


[g{Xs)-g{X^ji_^)]dBl 


s; \\g\\ocCH^,H^VH^{X- [Tk-l,Tk])VH^{B“ ■, [Tk-l,Tk]) 




(6.5) 


|5(A,n)AB5| ^ sup \g{X^)\Ge.T{rk - = 0^{d^-^). 

ue[o.T] 


Therefore (6.3) holds. 

Next we prove (6.4). Since 


A(")x.n=/ f{Xs)ds- "f{Xs)ds+ g{Xs)dB^- 


g{Xs)dBt 


and 


r g{Xs)dB^- r~\{XMB^ = r [g{Xs)-g{X,j^ J]dBf 

J-S-i J-S-2 

- r~' [5(A,) - 5(A.n J]dBf + 5(A.n_ J (aB^ - AB?._] 


then by mean value theorem with K equal to Lipschitz constant of / 

fiXs)ds- f{X,)ds 

'k-l '''^fe-2 I 

i^2Kd„ max sup |X, - X.n I ^®= ’ 0„ (d^+^"®). 

*^-1' ^ ^ 

We conclude that (6.4) also holds. 


s; \f{Xs) - /(X.n_j| ds + \f{x.^_,) - f{Xs)\ 


ds 


6.1.3 Application 

Results of subsection 6.1.2 imply proposition constituting the basis of the first example. 

Proposition 6.1. Assume that a model defined by (1.1) satisfies conditions of subsection 6.1.2. If in addition (2.2) 
holds, then Theorems 2.1-2.2 apply. 


6.2 Example 2 

Consider the Verhulst equation 


Xt=^+ [ {\Xs-Xl)dt + a [ XsdBf, e> 0 , tG[ 0 ,r]. 

Jo Jo 

f{x) = \x — is not Lipschitz. So we can’t use results of section 6.1. It was proved in [13] that this equation admits 
an explicit solution 

^ exp{At + crBi^} 


Xt = 


Moreover, from [13] one can get 

AAtn =Xtn 


1 + ^ /jj exp{As + aBp}ds 


A - + uABfn + O,. - 


, t€[0,T]. 


A(^>Wg +0. (;^5(iT7y)) 


Finally, explicit form of Xt implies existence of an a.s. finite r.v. c such that sup^gjg j,] ^ C a.s. Thus, one can use 

Theorem 2.1 if constant a is known and Theorem 2.2 in general situation. 
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